Electron-electron interactions in graphene bilayers 



S. Viola Kusminskiy^, D. K. Campbell^ and A. H. Castro Neto^ 



O 
O 
(N 



(N 



I 

O 

o 



> 

o 
cn 
p 

o 
oo 
O 



Department of Physics, Boston University, 590 Commonwealth Ave., Boston, MA 02215 



PACS 81.05.Uw - Carbon, diamond, graphite 

PACS 71 . lO.-w - Theories and models of many-electron systems 

PACS 51 .35.+a - Mechanical properties; compressibility 

Abstract. - We study the effect of electron-electron interactions in the quasiparticle dispersion 
of a graphene bilayer within the Hartree-Fock-Thomas-Fermi theory by using a four-bands model. 
We find that the electronic fluid can be described by a non-interacting like dispersion but with 
renormalized parameters. We compare our results with recent cyclotron resonance experiments in 
this system. 



Since graphene was isolated in 2004 [1] , it has attracted 
attention because of its possible application in all-carbon 
based electronic devices [2] and its connections to rela- 
tivistic field theory [3] . While there is strong theoretical [4] 
and experimental evidence [2,5] that single layer graphene 
(SLG) behaves as essentially a weakly interacting gas of 
two-dimensional (2D) Dirac particles, the situation in bi- 
layer graphene (BLG) is much less clear. Early theoretical 
studies have indicated that the SLG is much less prone 
towards magnetic states [6], while BLG can become mag- 
netic at low densities [7]. Moreover, while the electronic 
compressibility of SLG has essentially features of an in- 
sulator [5,8,9], the BLG compressibility is, unlike the 2D 
electron gas (2DEG) [10], non-monotonic and strongly de- 
pendent on electronic density [11]. It has also been argued 
that, unlike SLG, BLG should be unstable towards many- 
body states such as a pseudospin magnet [12], a Wigner 
crystal [13], and an excitonic superfluid [14]. It has been 
demonstrated that BLG is a tunable gap semiconductor 
by application of a transverse electric field [15, 16], leading 
to extra flexibility in dealing with its electronic proper- 
ties [17,18]. While electrons in BLG have a different topo- 
logical (Berry's) phase than electrons in SLG, as evident in 
integer quantum Hall effect measurements [19], the exper- 
imental evidence for electron-electron interaction effects in 
BLG has been elusive. Nevertheless, recent cyclotron res- 
onance experiments in bilayer graphene [20] have shown 
departures from the non-interacting bilayer model pro- 
posed by McCann and Falko [22]. These disagreements 
do not seem to be describable in terms of disorder effects 
alone [23]. The objective of our paper is to clarify these 
discrepancies. 



The SLG has a honeycomb lattice structure that leads 
to a Dirac-like electronic dispersion, £'(k) = ±c|k|, at 
the edges (the K and K' points) of the Brillouin zone. 
The electrons are described in terms of a 2D "relativistic" 
Dirac Hamiltonian with zero rest mass, where the velocity 
of light, c, is replaced by the Fermi-Dirac velocity, c. In 
the BLG (Bcrnal structure) the two graphene layers are ro- 
tated by a relative angle of 7r/3 that breaks the sublattice 
symmetry leading to 2 pairs of massive Dirac particles at 
the K (K') point. Nevertheless, the system remains metal- 
lic because 2 bands, belonging to different pairs, touch in 
a point. More explicitly, the non-interacting bands have 
the form: £;i(k) = -mS^ + £;(k), £:2(k) = mc^ - E{]s.), 



E{k) and £;4(k) = -mS^ - £;(k), where 



E{k) = ^(m52)2 + (5k)2. Hence, Ei{k) and E^ik) 
(i?2(k) and £^3(k)) describe a massive relativistic disper- 
sion with rest mass energy given by mc^. Rotations by 
other angles do not break the sublattice symmetry and 
hence do not lead to mass generation [24]. 

Our results suggest that BLG behaves as a liquid of 
Dirac quasiparticles with renormalized mass and velocity. 
The situation described here is unique when compared to 
standard non-relativistic Fermi liquids such as '^He [25] 
and ordinary metals [26] , or even to relativistic Fermi liq- 
uids such as quark matter in the core of neutron stars [27]. 
While the electrons in graphene are effectively "relativis- 
tic" , in the sense that they obey an effective Lorentz in- 
variance (only true at low energies) with the Dirac veloc- 
ity playing the role of velocity of light, on the other hand, 
from the point of view of an external observer, the whole 
graphene system is Galilean invariant and non-relativistic 
since the Dirac velocity is much smaller than the actual 
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speed of light. As a consequence, electron-electron inter- 
actions, just as in the case of relativistic [27] and non- 
relativistic [25] Fermi liquids, renormalize the quasiparti- 
cle mass, but unlike the relativistic and non-relativistic 
Fermi liquids, the "velocity of light" is also renormal- 
ized [28]. Moreover, in BLG as well as in SLG the ex- 
istence of negative energy bands and a pseudo-spin degree 
of freedom has to be included when constructing a Fermi 
liquid theory. This has been analyzed microscopically in 
references [29,30] and through a phenomenological model 
in Ref. [31] for SLG, however a treatment of BLG which 
takes into account its full hyperbolic four-bands structure 
is missing. 

The Coulomb interaction between electrons breaks the 
effective Lorentz invariance of the non-interacting prob- 
lem since it can be thought as instantaneous from the 
point of view of the electrons. In SLG, this violation of 
Lorentz invariance leads to the famous upward logarith- 
mic renormalization of the Fermi-Dirac velocity originally 
proposed in Ref. [32] . That effect is a result of the lack of 
screening in the SLG due to the vanishing of the density 
of states. The BLG, however, has a finite density of states 
at the Dirac point and hence screening plays an impor- 
tant role [33]. We show that, similarly to the 2DEG [26], 
the Hartree-Fock (HF) theory alone leads to an unphysi- 
cal logarithmic singularity at the Fermi surface indicating 
the importance of screening in this system. When screen- 
ing is accounted through the Thomas-Fermi (TF) theory, 
the log singularity is suppressed and, surprisingly, Lorentz 
invariance is recovered. We show that this result is due 
to the suppression of the intra-band transitions relative 
to the inter-band transitions. In order to test our find- 
ings, we study the problem of cyclotron resonance in this 
system and find good quantitative agreement with recent 
measurements [20]. 

We use a tight binding description of BLG in which 
only the in-plane, t « 3eV, and the out-of-planc, t± w 
0.37eV, nearest neighbor hopping parameters are consid- 
ered [35]. In this case we have c = 3ta/2 (a = 1.42 A) 
and mc^ = t±/2. From now on we choose units such that 
h ~ 1 = c. The hyperbolic shape of the non-interacting 
dispersion introduces an intrinsic energy scale in the prob- 
lem, m. In the "non-relativistic" (NR) limit, k ^ m, one 
can replace the four hyperbolic bands by two parabolic 
bands E±{k) = ±fc^/2m [22] . In this approximation, the 
usual NR dispersion of the 2DEG is recovered, but with 
allowed negative energy values. In the "ultra-relativistic" 
(UR) limit, » m, E{k) ~ fc, one obtains the SLG dis- 
persion. The crossover energy scale from non-relativistic 
to ultra-relativistic (NR-UR) is given by m. We notice 
that the effective non-relativistic low energy approxima- 
tion fails when treating the interacting problem because 
the Coulomb energy associated with electron-electron in- 
teractions is of the order of the inter-band transitions [11]. 

The electronic interactions are included by adding to 
the non-interacting energy Eo = I/tt^^-/ Ei{p)dp an 
exchange term which can be written as (energies are given 
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Fig. 1: Quasiparticle dispersion, ei(q) (in eV), as a function of 
the momentum q (in units of A). Inset: intra (i3+) and inter 
(D-) band contributions. Full: /3 = (HF); Dashed: /? = 1; 
Dash/dot:/3 = 2; Dash/dot-dot: /3 = 5; Dotted line: non- 
interacting. Curves have been shifted for comparison [kp / K = 
0.03 and g = 0.5). 

per unit area, and the spin and valley degeneracy factor 
of 4 is accounted): 

E,,=-^l xS(q,p)x"^(p,q)n.(9)^j(p)K.(q-p) (1) 

where a=itl correspond to the symmetric/antisymmetric 
representations of the Coulomb interaction: 

Kt(k) = 27re2(l ± exp{-fcd})/[2e(fc + /3m)] , (2) 

ni{q) is the occupation number of band i, and X'^ (q, p) are 
overlap matrices which contain information of the change 
of basis [7]. Screening is taken into account through 
the TF approximation by introducing a screening length 
in (2) that is proportional to the density of states. In 
(2) /3 is the parameter that controls the value of the 
TF screening length, the HF theory is obtained by tak- 
ing f3 ~ Within the Random Phase Approximation 
(RPA), Prpa = 4.9 (1 + E'f/™), being g /fic the 

dimcnsionlcss coupling constant, and Ep the Fermi en- 
ergy. For experimentally realized densities (jig « 10^^ - 
10^3 cm-2) and g = 0.5 (e = 3.9 for SiOa), /3 « 1 - 5. 
The energy of a quasiparticle in the i*'' band is given by 
e,{q) = (5i;/(5n,(q)|„^=„o, where (5n,(q) = n,;(q) - n°(q), 
being n^(k) the occupation number of the non-interacting 
system. E[5ni\ is the total energy E = Eq + Eex- We 
can therefore write {q) = Ei [q] + Ai?i [q) with Ai?.; {q) ~ 
-'^IqT,a,j Xg(q,P)x"i(P,q)"-°(p)K«(q-p) the correction 
to the non-interacting band Ei(q). 

We consider the case of electron doping such as that the 
chemical potential does not reach the uppermost band, 
which is usually the experimentally realized situation [36] . 
Therefore, our results are valid for Fermi energies up 
to \/2ix, which corresponds to densities smaller than 
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Fig. 2: v*{q) (in units of c) as a function of q (in units of the 
cutoff A) for kp/A = 0.03. Dotted line; non-interacting. Solid 
line: HF. Dashed line: TF with l3 = 1; Dash/dot: P = 2; 
Dash/dot-dot: (3 — 5. Inset: zoom-in near the Fermi vector 
and the logarithmic fit to the divergence (dashed line) for HF. 

He ~ 10^"^ cm"^. We look then at the correction to the first 
band, which we write as AEi{q) = D+{q, kp) + D-{q, A) 
to distinguish intra-band (£*+) from inter-band (-D-) con- 
tributions. The expressions for D± can be easily de- 
rived from AEi{q). kp is the Fermi wave vector and 
A a cut-off of the order of the inverse lattice spacing 
(A « lA"^ « 7eV). 

Fig. 1 shows the quasiparticle band within the HF the- 
ory (solid line) for a typical value of the Fermi vector. Fig. 
1 (a) depicts the correction due to the intra-band transi- 
tions, £)_(.. Its behavior, as expected, is qualitatively very 
similar to that of a 2DEG [26]. In particular, the inflection 
point seen at q ~ kp is due to the special role of kp which 
separates a domain with an avoidable singularity q < kp 
from a singularity free domain for q > kp. While Dj^- 
diminishcs with kp, the correction due to is indepen- 
dent of it. The latter is shown in Fig. 1 (b). As it can be 
seen from the figure, for typical electronic densities, the 
correction due to inter-band interactions is roughly two 
orders of magnitude bigger than that of the intra-band. 
Notice, from Fig. 1, that the quasiparticle dispersion how- 
ever inherits the inflection point from at q ~ kp. The 
renormalized band velocity is given by v*{q) = \de/dq\, 
which is plotted in Fig. 2. Due to the sharp inflection 
point in AEi at kp, the effective quasiparticle Fermi veloc- 
ity, v*{kp) presents an unphysical logarithmic divergence: 
v*{k ~ kp) ^ —4g/Tr\og{\k — kp\/A), as it occurs for 
the 2DEG. For small momentum nevertheless, q/kp ^ 1, 
the renormalized dispersion can be shown to be parabolic: 
ei{q) w q^/{2m), with m^^ = m^^ + g{5/kp - l)/2. 

As mentioned earlier, the divergence of the Fermi ve- 
locity is an unpleasant feature of the HF approximation 
which indicates the necessity of introducing screening in 
the problem. The renormalized band ei (g) is shown in Fig. 
1 for different values of (3. We sec that the introduction of 
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Fig. 3: Left: m* (in units of the bare mass m) as a function of 
electron density (in units of electrons per cm~'^). Right: c* 
(in units of the bare velocity c) as a function of rie, for different 
P (3 = 0.5). 

screening eliminates the inflection point at kp. This can 
be seen clearly in Fig. 2, where it is shown that the diver- 
gence in the quasiparticle velocity disappears for finite p. 
The most striking feature of our calculations is that the 
quasiparticle dispersion can be fitted by a non-interacting 
like dispersion: 

ei{k)^eo + ^{v*pky + {m*c*^)\ (3) 

where eo is a constant, and m* and c* arc the quasipar- 
ticle mass and renormalized "light" velocity, respectively. 
We find that this result is valid to high accuracy for a 
large region of energy and momenta due to the fact that 
the inter-band transitions largely dominate over the intra- 
band ones [37]. 

The results for m* /m and c*/c are shown in Fig. 3 as a 
function of electronic density. He = kp/ir, for different val- 
ues of the screening strength /3. Note that, for fixed (3, c* /c 
increases monotonically with density, whereas m* /m has 
a minimum at a finite rig. While m* /m is renormalized 
to smaller values, c*/c is renormalized to larger values. 
This has interesting consequences for the NR-UR crossover 
mentioned earlier. The crossover energy for the non- 
interacting problem is given by = mc^. Analogously, 
for the interacting result we can define the crossover en- 
ergy as El = m*{c.*Y. This quantity is plotted in Fig. 4 
(b), which shows that E* > Ec for all the values of the 
parameters (this is also true if we vary the coupling con- 
stant 0.1 < 5 < 2). However, the relevant parameter to 
compare with experiment is the crossover electronic den- 
sity. Tie = q^/T^ = (TOc)^/7r, that is, the density at which 
the NR-UR crossover takes place. Fig. 4 (a) shows the 
renormalized value of this quantity, n* = (m*c*)^/7r, in 
units of the non- interacting value ric. Indeed, it is seen 
that n* < Uc always, even though the renormalized quasi- 
particles' energy is higher. 

Let us now consider the problem in the presence of a 
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Fig. 4: Left: NR-UR crossover density n* (in units of the bare 
crossover density ric) as a function of electronic density (in 
units of electrons per cm~ ). Right: (in units of the non- 
interacting crossover energy Ec) as a function of rie, for differ- 
ent P {g = 0.5). 

transverse magnetic field B. For the non- interacting prob- 
lem, the Landau levels are given by (restoring units) [38]: 
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l + 16rVl6r2 
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(4) 



UJc 
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c^/leB jc is the cy- 
One can clearly see 



where n is a positive integer 
clotron frequency, and r ~ mfv^^juj 
that this problem has the NR-UR crossover as a func- 
tion of B discussed earlier. At low fields, r ^ 1, we find 
« ±[w^/(2toc^)] \Jn{n ■\- 1), and the Landau level en- 
ergy is proportional to B as in the NR problem [22]; at 
high fields, r <C 1, one finds k, ±tJc \fn and, as in the 
UR case, we find the Landau level energy proportional to 

4b. 

Just as in the case of a Fermi liquid, here the quasi- 
particles carry electric charge e and couple to a magnetic 
field via minimal coupling. Note that here, however, the 
cyclotron mass is not protected by Kohn's theorem since 
the dispersion is not parabolic [39,40]. Hence, the Landau 
level spectrum is the same as the non-interacting problem, 
eq. (4), with the bare parameters, m and c, replaced by 
renormalized ones, m* and c*, respectively. In Fig. 5 we 
show the data from cyclotron resonance experiments [20] 
for inter-Landau level transitions for different filling fac- 
tors V = ne/ng> (n$ is the density of flux quanta through 
the system) together with our results for g = 0.5, /3 = 4, 
5 = 0.76x10^ m/s (c* = 1.2-1.3x10*^ m/s) andt_L « 0.33 
eV {t*j_ = 0.44 — 0.45 eV). These values were obtained by 
fitting our renormalized theory to the experimental data, 
taking /3, g, and the bare m and c as free parameters, 
and are in agreement with recent infrared spectroscopy 
data [21]. The variation of c* and m* with electronic den- 
sity was taken into account, however we took the screening 
strength /3 as fixed, since our treatment is not self consis- 
tent. Nevertheless, as a double check, the value obtained 
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Fig. 5: Landau level transition energies for different filling fac- 
tors 1/ = 4 (n: -> 1), ^ = 8 (n: 1 ^ 2), u ^ 12 (n: 2^3) 
and = 12 (n: 2 —r 3) as a function of B (in Tesla). The 
experimental data (dots) was taken from Ref. [20]. The solid 
line is the theoretical value for g = 0.5, and /3 = 4 for bare 
values c ^ 0.76 x 10*^ m/s and t± x 0.33 eV. The dashed line 
shows the fits obtained in Ref. [20] for t± = 0.35 eV. 



for g is the expected for BLG on Si02 and j3 falls within 
the expected range for such densities. One can see that our 
results are in fair quantitative agreement with the exper- 
imental data, giving support to the idea that this system 
can be described by a Dirac liquid of quasiparticles with 
a dispersion given by (3). There is a small electron- hole 
asymmetry due to inter-band interactions (and therefore 
independent of density) , which results in smaller values of 
t*j_ (~ 10%) and c* (~ 5%) for hole doping. This difference 
however is not enough to explain the asymmetry observed 
in [20]. As it was mentioned above, the screening strength 
(3 was fitted to a constant value for all four plots in Fig. 
5, and it should be taken as the best average (3 that gives 
a reasonable good fit for all the data range. For a better 
agreement with the data it would be probably necessary 
to include selfconsistently the dependence of f3 with the 
electronic density but this goes beyond the scope of the 
present paper. 

We have studied the effect of electron-electron interac- 
tions on the electronic properties of a graphene bilayer 
within the Hatree-Fock-Thomas-Fermi theory by taking 
into account the full four-bands model of BLG. We have 
shown that the quasiparticles can be described by a non- 
interacting, Lorentz like dispersion with renormalized pa- 
rameters which depend on the electronic density. The fact 
that the Lorentz invariance of the dispersion is recovered 
for a large range of energies is an unexpected result since 
there is no evident symmetry behind it. It is important to 
note that this result is due to the dominance of inter-band 
transitions, which are missing in the usual Fermi liquid 
picture. Since this contribution is independent of elec- 
tronic density, the accuracy of the effective description 
increases with decreasing density. Further investigation 
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is needed to determine if corrections beyond HF lead to 
deviations from the Lorentz dispersion. Furthermore, we 
have tested om" calculations by comparing our results with 
recent cyclotron resonance experiments and found quan- 
titative agreement between theory and experiment. Our 
results are also in agreement with recent ab-initio calcula- 
tions [41]. 

The authors wish to thank A. K. Geim, Z. Q. Li, A.H. 
MacDonald and Johan Nilsson for their valuable com- 
ments on the manuscript. 
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